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ABSTRACT 


This  report  presents  the  results  of  simultaneous  HF  Doppler  and  UHF  in¬ 
coherent  backscatter  measurements  of  the  travelling  ionospheric  disturbances 
observed  at  Arecibo  (P.R.)  during  the  summer  of  1967.  One  of  the  prime  objec¬ 
tives  in  conducting  this  series  of  experiments  involving  3  spaced  phase  path 
sounders  and  the  backscatter  radar  is  to  examine  the  gravity  wave  interpretation 
of  the  large  scale  travelling  disturbances.  The  analysis  of  the  observations  is 
based  upon  three  distinct  techniques:  (1)  Cross-correlation  analysis  of  the  HF' 
Doppler  records  for  an  accurate  estimate  of  the  wave  velocity,  (2)  HF  ray  tracing 
for  demonstrating  the  self-consistency  of  the  HF  and  the  incoherent  backscatter 
data,  and  (3)  the  coupling  of  acoustic  gravity  waves  to  the  ionospheric  plasma. 

It  is  concluded  from  the  results  of  the  analysis  that  the  large  scale  travelling 
ionospheric  disturbances  are  truly  a  manifestation  of  the  acoustic-gravity  waves. 
Finally,  some  aspects  of  the  analysis  which  will  be  refined  in  the  next  phase  of 
the  study  are  also  outlined. 
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SECTION  I.  INTRODUCTION 


This  report  is  concerned  with  the  question,  "Are  travelling  ionospheric 
disturbances  acoustic  -  gravity  waves  or  are  they  not?"  The  development  of 
acoustic  -  gravity  wave  theory  in  recent  years  has  made  it  clear  that  if  meaning¬ 
ful  comparisons  between  theory  and  observation  are  to  be  made,  the  velocity, 
period,  and  vertical  structure  of  the  disturbances  must  be  measured  simultaneous¬ 
ly.  ^  HF  radar  measurements  over  the  past  20  years  have  shown  that  wave 
periods  from  a  few  minutes  to  a  few  hours  and  propagation  velocities  from  less 
than  100  to  more  than  1000  meters  per  second  can  be  expected.  ^  All  these  can 
be  explained  as  acoustic  gravity  waves  if  suitable  assumptions  about  the  vertical 
structure  are  made.  Similarly,  the  period  and  vertical  structure  of  a  few 

(3) 

travelling  disturbances  have  been  measured  by  means  of  incoherent  backscatter 

and  these  observations  are  consistent  with  the  predictions  of  acoustic  -  gravity 

(4) 

wave  theory  so  long  as  suitable  velocities  are  assumed. 

In  this  report  a  series  of  experiments  is  described  in  which  simultaneous 
HF  radar  and  incoherent  backscatter  observations  of  travelling  ionospheric  dis¬ 
turbances  are  made  and  the  need  for  assuming  various  wave  parameters  elimi¬ 
nated. 

The  report  begins  with  a  discussion  of  the  factors  which  influenced  the  ex¬ 
periment  design,  followed  by  a  brief  description  of  the  equipment  involved.  The 
observations  are  then  summarized  and  selected  pieces  of  experimental  data  are 
shown.  An  outstanding  example  of  a  travelling  ionospheric  disturbance  is  selected 
and  analyzed  in  detail.  Finally,  conclusions  are  drawn  and  the  work  planned  for 
the  final  phase  of  this  study  is  outlined. 
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SECTION  II.  EXPERIMENT  DESCRIPTION 


A.  Station  Geometry 

The  experiment  w?.s  designed  to  study  what  are  commonly  called  "large 
scale"  travelling  ionospheric  disturbances.  Typically  these  events  have  time 
scales  of  10-30  minutes  and  propagation  velocities  on  the  order  of  100  meters 
per  second.  They  occur  at  the  rate  of  about  ie  per  day  and  cause  electron 
density  perturbations  of  perhaps  10  percent.  Shorter  period  disturbances  (a  few 
minutes)  occur  more  frequently  but  cause  much  smaller  density  perturbations  and 
therefore  are  difficult  to  study  with  incoherent  backscatter.  Longer  period  dis¬ 
turbances  (a  few  hours)  cause  much  larger  density  perturbations  but  occur  in¬ 
frequently,  perhaps  once  a  month.  In  view  of  the  competition  for  observing  time 
at  incoherent  backscatter  facilities,  the  very  small  probability  of  occurrence 
for  these  long  period  events  makes  their  study  impractical. 

The  horizontal  velocity  of  propagation  is  measured  by  spacing  three  identi¬ 
cal  HF  sounders  on  the  island  of  Puerto  Rico  and  measuring  time  delay  between 
the  appearance  of  the  disturbance  at  each  station.  The  spacing  between  sounders 
is  a  compromise.  The  spacing  must  be  large  enough  so  that  significant  time 
delays  will  be  observed  but  close  enough  so  that  the  waveform  of  the  disturbance 
will  not  change  appreciably  from  station  to  station.  Experimenters  in  the  past 
have  found  a  spacii._,  of  30-50  km  near  optimum  for  disturbances  of  the  time  scale 

(5' 

being  sought  here.  '  Convenient  locations  were  found  at  the  Arecibo  Ionospheric 
Observatory  where  the  incoherent  backscatter  observations  were  made,  near 
Mayaguez  where  Prof.  Braulio  Dueno  mai.  tains  a  radio  field  station,  and  at  the 
Ramey  Air  force  Base.  The  geometry  of  network  is  shown  in  Figure  1. 

B.  Equipment  Characteristics 
1 .  HF  Rounder  Network 

Each  of  the  three  HF  sounder  installations  consisted  of  two  phase-coherent 
pulse  radars  which  could  be  operated  throughout  the  3-30  MHz  band.  In  practice 
the  radar  frequencies  were  chosen  to  lie  below  the  vertical  incidence  critical  fre¬ 
quency  and  all  three  sounders  were  operated  on  the  same  pair  of  fiequencies.  A 
pulse  width  of  200  microseconds,  a  pulse  repetition  frequency  20  pps,  and  a  peak 
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pulse  power  of  20  kw  were  used  a*  each  site.  The  transmitted  frequencies,  the 
receiver  local  oscillator  frequencies,  and  the  timing  circuitry  were  controlled 
from  a  stable  master  oscillator.  The  master  oscillators  at  all  sites  were  syn¬ 
chronized  to  avoid  interference  between  stations.  The  receiver  output  was  beat 
to  a  5  cps  IF  £•  recorded  on  analog  mangetic  tape.  The  tapes  were  mailed  to 
the  Raytheon  laboratories  at  Burlington,  Massachusetts  for  spectrum  analysis. 

2.  Incoherent  Backscatter  Facility 

The  incoherent  backscatter  observations  were  made  using  the  '*30  MHz  rada 
at  the  Arecibo  Ionospheric  Observatory.  Sixty  microsecond  pulses  were  transmit 
ted  with  a  peak  power  of  2.  5  megawatts.  The  pulse  repetition  frequency  was  90 
pulses  per  second.  The  antenna  beam  was  vertically  oriented.  At  F- region 
heights  the  beam  is  about  a  kilometer  in  diameter  so  that  the  scattering  volume 
in  space  was  a  cylinder  of  9  kilometers  long  and  about  one  kilometer  wide.  Both 
the  ion  and  electron  components  of  the  backscatte red  power  spectrum  were  re¬ 
corded  on  digital  tape. 
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SECTION  III.  OBSERVATIONS 

A.  Data  Collection 

The  network  of  three  ph ise-oath  sounders  was  available  for  the  months  of 
May,  June,  and  July  1967.  During  this  time,  observations  were  made  whenever 
time  on  the  430  MHz  radar  could  be  obtained  during  daylight  hour?  (0700-1900 
local  time).  Daytime  operation  was  required  because  the  electron  component  of 
the  back- scattered  power  spectrum  (the  "plasma  line")  is  only  detectable  when  it 
is  enhanced  by  the  daytime  flux  of  photoelectrons.  In  all,  263  hours  of  observing 
time  were  scheduled  and  of  this,  all  equipment  worked  satisfactorily  during  185 
hours. 

B.  Data  Reduction 

All  usable  data  have  been  reduced  and  examined.  The  incoherent  backseat- 
ter  observation?  were  reduced  at  AIO  to  contour  plots  of  constant  plasma  fre¬ 
quency  versus  height  and  time.  The  HF  observations  on  both  frequencies  and 
from  each  station  were  reduced  at  Raytheon  to  spectrograms  of  Doppler  shift 
versus  time  for  various  ran~e  gate  positions.  Events  of  special  interest  were 
recycled  through  the  data  reduction  process  one  or  more  times  to  optimize  the 
presentation. 

C.  Results 

1 .  Summary 

Travelling  ionospheric  disturbances  are  found  to  be  a  common  feature  of  the 
daytime  ionosphere.  Roughly  half  of  the  HF  Doppler  data  contain  anoma.ies  which 
can  be  followed  from  station  to  station  with  little  change  in  form.  The  statistics 
of  the  observations  will  be  discussed  in  some  detail  in  the  final  report  where  cor¬ 
relations  with  possible  geophysical  sources  will  be  sought.  For  the  purposes  of 
this  report,  the  data  as  a  whole  is  summarized  simply  by  tabulating  the  speed  and 
direction  of  all  prominent  events.  The  velocities  have  been  computed  from  hand- 
scaled  time  delays  between  stations.  The  results  are  shown  in  Table  1.  The  com¬ 
puted  velocities  are  quite  comparable  to  those  found  in  earlier  studies.  ^  It  is 
concluded  that  the  disturbances  observed  in  this  study  are  of  the  same  type  as 
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Table  1.  Horizontal  Velocities  of  the  Ionospheric  Disturbances 


Date 

Time 

(U.T.) 

Speed 

(M/S) 

Bearing 

(Deg.  West  of  North) 

5-07-67 

Hr.  Min. 

13  :  14 

66.61 

4.  36 

5-1 ?-67 

1  i  :  50 

225. 52 

55.  16 

51  4—67 

11 

35 

218. 32 

17.  60 

5-18-67 

21 

40 

148.  10 

231. 21 

5-20-67 

13 

30 

169. 47 

83.  40 

5-21-67 

12 

25 

154. 30 

76.  20 

5-21-67 

21  :  20 

218.  32 

197. 60 

5-25- 57 

22  :  45 

679. 50 

55.  20 

5-27-67 

21  : 

05 

680. 00 

55.  00 

5-30-67 

22  : 

:  55 

428. 30 

141. 16 

6-01-67 

i  3  ;  48 

197. 74 

60.  44 

6-C3-67 

22  : 

:  38 

85.73 

102. 45 

6-08-67 

13  : 

:  45 

374.  74 

159. 00 

6-G8-67 

22  :  49 

124. 31 

191. 20 

6-14-67 

12 

00 

112. 70 

124. 50 

6- -5-67 

12 

39 

172. 41 

129. 50 

6-15-67 

13 

35 

124. 46 

41.76 

6-18-67 

22 

48 

177.  50 

101. 00 

6-19-67 

22 

54 

1 79. 04 

149. 18 

6-24-67 

23 

00 

226.07 

177. 84 

6-27-67 

1 

10 

277. 58 

170.  * 

6-27-67 

22 

10 

122. 75 

169.  10 

6-23-67 

16  :  55 

138. 79 

350. 44 

6-26-67 

20 

52 

180. 50 

199. 14 

6-30-67 

13 

30 

107. 20 

126. 58 

7-07-67 

11 

57 

130. 82 

"8.  16 

7-08-67 

11 

12 

141. 43 

113. 99 

7-12-67 

13 

38 

159. 53 

172. 17 

7-13-67 

12 

32 

275. 81 

15.  24 

7-15-67 

13 

25 

161. 41 

186.47 

/  - 16-67 

11 

50 

132. 24 

37.  24 
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Table  1  Cont. 


Date 

Time 

(U.T.) 

Speed 

(M/S) 

Bearing 

(Deg.  West  of  Nor 

7-21-67 

Hr.  Min. 

22  :  22 

186. 46 

191.20 

7-22-67 

13  :  02 

191.69 

34.  25 

7-29-67 

14  :  15 

155. 55 

107.  93 

7-31-67 

19  :  33 

201. 42 

159. 00 

8-01-67 

22  ;  lo 

187. 28 

174. 52 

th) 
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those  studied  extensively  in  the  past.  A  t.  st  of  the  acoustic-gravity  wave  inter¬ 
pretation,  based  on  the  data  iecently  collected  is  therefore  believed  applicable  to 
earlier  observations  as  well. 

2.  Examples 

Three  cases  were  observed  of  what  may  be  called  "classical"  travelling 
ionospheric  disturbances.  These  are  classical  in  the  sense  that  when  HF  experi¬ 
menters  wish  to  give  an  example  of  their  data,  a  disturbance  of  this  type  is  usual- 
ly  shown.  ’  These  events  are  characterized  by  "switch  backs"  on  plots  of 
Doppler  shift  versus  time  (or  group  path  versus  time  if  this  parameter  is  analyzed). 
While  these  events  are  more  spectacular  than  most,  it  is  clear  that  they  differ 
from  disturbances  which  produce  single-valued  Doppler  signatures  cnly  in  fhat 
the  horizontal  gradients  of  electron  density  they  cause  are  slightly  greater.  This 
usually  means  that  the  wave  is  high  in  amplitude  and  therefore  will  be  easily 
recognized  with  incoherent  backscatter.  Multiple  valued  disturbances  have  there¬ 
fore  been  given  special  attention  in  this  report. 

The  first  example  is  shown  in  Figure  2.  There  are  two  major  disturbances 
here,  one  a*  1410-1425  and  one  at  1430-1445.  The  Doppler  3hift  is  multiple 
valued  during  these  periods,  indicating  that  the  contours  of  constant  electron 
density  in  the  ionosphere  are  distorted  enough  so  that  multiple  reflections  take 
place.  Time  shifts  of  a  few  minutes  are  observed  between  stations.  The  corres¬ 
ponding  incoherent  backscatter  data  is  shown  in  Figure  3.  The  contours  of  con¬ 
stant  plasma  frequency  are  derived  directly  from  the  ion  component  of  backncat- 
tered  power  and  normalized  to  the  critical  frequency  at  the  layer  peak  as  read 
from  ionograms  taken  simultaneously.  The  ionosphere  is  clearly  affected  by 
waves  over  a  considerable  height  range  including  the  level  of  HF  reflection  (about 

200  kilometers  for  5.228  MHz).  The  phase  fronts  oi  the  wave.,  are  tilted  forward 

(3) 

as  expected  from  earlier  incoherent  backscatter  results. 

The  second  example  is  shown  in  Figure  4.  The  time  scale,  Doppler  shift 
structure,  and  time  delays  between  stations  are  similar  to  the  event  just  discussed. 
Unfortunately,  the  430  MHr  r  darat  AIO  was  not  operable  during  this  event  and 
the  incoherent  backscatter  data  is  not  available. 
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Figure  3.  Contours  of  Constant  Plasma  Frequencies  from  Incoherent 
Backscatter  Measurements  at  430  MHz  on  June  12,  iv67 
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Figure  4.  Doppler  Shift  Recordings  on  June  14,  1967  at  an  Operating 

Frequency  of  6.843  MHz 


The  third  example  is  shown  in  Figure  5.  The  classical  events  occur  about 
1130  and  1125.  The  data  from  1130  to  1400  is  show’n  because  it  corresponds  to  a 
major  wave  event  on  the  incoherent  backscatter  data.  The  backscatter  data  is 
shown  in  Figure  6.  The  diamonds  and  connecting  lines  between  250  and  300  kilo¬ 
meters  indicate  the  height  of  the  F- region  peak.  The  "islands"  at  the  peak  which 
form  at  about  1140,  1225,  and  1305  reprerent  condensations  caused  by  a  train  of 
travelling  waves.  The  sloping  phase  fronts  of  the  wave  crests  are  evident. 

During  much  of  this  disturbance,  the  BF  ray  paths  were  severely  distorted  and 
reliable  measurements  of  wave  velocity  can  not  be  made  using  hand  scaled  delay 
times  between  stations.  It  is  expected,  however,  that  the  cross- correlation 
analysis  which  would  be  extended  to  this  data  during  the  next  phase  of  the  study 
will  provide  meaningful  measurements  of  the  velocity  throughout  this  event. 
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Figure  5.  Doppler  Shift  Recordings  on  July  8,  1967  at  an  Operating 

Frequency  of  5.  328  MHz 
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Figure  6.  Contours  of  Constant  Plasma  Frequencies  from  the 
Backscatter  Measurements  on  July  8,  1967 


SECTION  IV.  ANALYSIS 

The  event  of  12  June  1967  is  representative  of  the  J  classical  travelling 
ionospheric  disturbances  observed  and  will  be  used  for  a  comparison  with  acoustic- 
gravity  wave  theory.  The  analysis  is  based  upon  three  distinct  techniques:  Cross 
correlation  for  an  objective  measurement  of  the  wave  velocity,  ray  tracing  for 
demonstrating  the  self-consistency  of  the  HE  and  the  incoherent  backscatter  data, 
and  the  coupling  of  acoustic  gravity  waves  to  the  ionospheric  plasma.  The  tech¬ 
niques  themselves  are  discussed  in  the  appendix  and  referred  to  in  the  analysis 
when  used. 

The  HF  and  incoherent  backscatter  data  for  June  12  have  been  shown  in 
Figures  2  and  3  respectively.  The  first  step  in  the  analysis  is  to  derive  the  hori¬ 
zontal  velocity  of  the  wave  by  cross  correlation  of  the  HF  signatures.  The  HF 
signatures  have  been  dig  tized  by  hand  to  provide  a  one-bit  representation  for  the 
computer.  Figure  7  shows  the  digitized  signals  at  the  three  stations.  The  cross¬ 
correlation  program  is  described  in  Appendix  A.  Figure  8  shows  the  result  for 
the  20  minute  window  containing  the  first  and  largest  signature.  The  derived 
speed  is  112  m/s.  The  direction  of  travel  is  from  75*  west  of  north. 

The  next  step  in  the  analysis  is  to  use  the  derived  speed  to  convert  the  time 

axis  of  the  incoherent  backscatter  data  (Figure  ?)  to  horizontal  distance.  This 

defines  a  two-dimensional  distorted  ionosphere  which  moves  overhead  at  a  known 

speed.  The  assumption  of  two-dimensionality  is  equivalent  to  assuming  that  the 

disturbance  is  a  plane  wave  and  past  studies  have  shown  this  to  be  at  least  ap- 

(2) 

proximately  true. 

At  this  point  a  check  on  the  self-consistency  of  the  experimental  data  has 
been  made  by  doing  ray  tracing  through  the  distorted  ionosphere  to  synthesize  the 
HF  signature  observed  at  Arecibo.  The  ray  tracing  program  is  described  in 
Appendix  B.  The  result  is  shown  in  Figure  9.  The  match  is  felt  to  be  surprisingly 
good  considering  the  many  idealizations  made  to  this  point  and  it  is  concluded  that 
the  experimental  data  is  self-consistent. 
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TIME-GMT 


Figure  7.  One-bit  Digitized  Version  of  the  Doppler  Shift  Signatures 

of  June  1?.,  1%7 


CORRELATION  COEFFICIENT 


14  08  30-  14  ■■  28  30  UT 


Figure  8.  Cross  Correlograms  for  the  3  Pairs  of  Stations  for  a  Selected 
20  Minute  Window  Centered  on  the  Main  Signature  of  the  Digitized 

Doppler  Record 
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Figure  10  compares  the  reflection  'evel  for  the  HF  radar  (ordinary  ray)  as 
derived  from  incoherent  backscat^er  with  the  HF  signatures  observed  at  Arecibo. 
Clearly,  tf  e  signature  is  associated  with  the  concave  ripples  in  the  contours.  The 
wave  period  is  about  20  minutes  and  the  wave  length  about  1^6  kilometers.  Figure 
11  shows  the  level  of  reflection  with  the  appropriate  distance  scale  added.  Also 
shown  is  a  line  of  constant  phase  predicted  by  acoustic-gravity  wave  theory.  The 
equation  shown  on  this  figure  is  the  dispersion  relation  for  such  waves  where  K 

Z 

is  the  (complex)  vertical  wave  number,  co^  is  the  angular  brunt  frequency  of  the 
neutral  atmosphere,  u>  is  the  angular  wave  frequency,  is  the  (real)  horizontal 
wave  number,  H  Is  the  neutral  scale  height,  and  is  the  kinematic  viscosity  of 
the  neutral  gas.  This  equation  neglects  phase  shifts  due  to  production  and  recom¬ 
bination  but  is  correct  to  first  order.  The  derivation  including  these  processes  is 
given  in  Appendix  C.  Theory  predicts  that  at  all  altitudes  the  trough  of  the  wave 
should  fall  along  the  dashed  line. 

Contours  of  constant  plasma  frequency  over  the  entire  height  range  are 
shown  in  Figure  12.  Agreement  between  the  position  of  the  wave  trough  and  the 
dashed  line  is  seen  to  be  very  good.  Qualitatively,  the  wave  amplitude  should 
grow  exponentially  with  altituue  antil  viscous  damping  becomes  important.  As 
damping  sets  in,  the  lines  of  constant  phase  should  swing  towards  the  vertical  and 
the  wave  amplitude  should  drrp.  The  swing  to  the  vertical  is  observed  to  occur 
in  the  range  200-250  kilometers  altitude  and  the  amplitude  of  the  disturbance  is 
also  observed  to  drop  in  the  same  interval.  The  agreement  between  acoustic- 
gravity  wave  theory  and  the  observed  travelling  disturbance  structure  is  taken  as 
strong  evidence  for  the  correctness  of  this  interpretation. 
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430  MHz 


Comparison  of  Variations  in  5.228  MHz  Reflection  Level 
he  HF  Signature  at  the  Same  Frequency  at  Arecibo 


12  JUNE  1967  430  MHz 
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HORIZONTAL  DISTANCE- (K  M  ) 

Figure  11.  Predicted  Variation  of  the  Phase  of  an  Acoustic  Gravity 

Wave  with  Altitude 


SECTION  V.  CONCLUSIONS 


It  is  concluded  that  large  scale  travelling  ionospheric  disturbances  are 
acoustic-gravity  waves.  No  significant  discrepancies  1  'tween  theory  and  observa¬ 
tion  have  been  found.  In  the  final  phase  of  this  study,  some  aspects  of  the  analysis 
will  be  refined,  such  as  including  loss  and  production  terms  when  coupling  the 
neutral  wav  to  the  ionospheric  plasma.  It  is  clear,  however,  that  while  such 
refinements  may  point  up  areas  which  need  more  theoretical  work,  the  acoustic- 
gravity  wave  interpretation  itself  is  on  firm  ground. 

Future  plans  for  the  remaining  portion  of  this  study  include  three  tasks. 
First,  the  comparison  between  theory  and  observation  will  be  refined  as  described 
above.  Second,  the  assumption  (usually  made  in  HF  studies)  that  all  important 
Doppler  shift  anomalies  are  introduced  at  or  near  the  reflection  height  will  be 
examined.  Third,  a  search  fcr  the  wave  sources  will  be  made  by  looking  for 
correlations  between  the  observed  events  and  other  geophysical  disturbances 
occurring  at  the  proper  time  and  place. 
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SECTION  VI.  APPENDIXES 


A.  Cross  Correlation 

Horizontal  velocity  of  a  travelling  ionospheric  disturbance  is  determined 
by  measuring  time  delays  in  the  occurrence  of  similar  complex  signatures  on  the 
Doppler  shift  records  obtained  at  the  three  spaced  sites  as  stated  in  Section  III  (c). 
The  determination  of  the  time  lags  by  lining  up  the  similar  signatures  on  the  3 
Doppler  shift  records  involves  a  certain  degree  of  subjective  error  especially 
when  the  signatures  extend  considerably  in  time  and  are  not  of  identical  configura¬ 
tion  on  the  3  tracks.  A  statistical  approach  involving  the  calculation  of  cross 
correlation  functions  between  each  pair  of  the  Doppler  records  is  desirable  in 
situations  like  this  to  extract  the  best  estimates  of  the  time  delays  and  thereby 
the  required  velocity  information.  The  cross  correlation  analysis  is,  in  fact, 
wii  :ly  used  for  studying  the  structure  and  movement  of  the  small  scale  irregu¬ 
larities  om  the  observed  3  station  amplitude  fading  patterns  of  a  reflected  radio 
signal  '.  The  Doppler  data  collected  on  June  12,  1967,  which  exhibited  a 
distinct  loop  type  complexity  at  the  3  stations,  is  subjected  to  the  correlation 
analysis  as  described  below.  Unlike  the  amplitude  patterns,  the  Doppler  shift  of 
a  received  signal  can  have  multiple  values  at  any  given  instant  of  time  and  this 
in  priciple  leads  to  a  correlation  function  in  the  two  dimensions  of  time  and  fre¬ 
quency.  In  the  present  analysis,  however,  the  correlation  functions  are  calculated 
in  the  single  dimension  of  time  only  for  the  case  of  zero  relative  shift  in  the 
frequency  (i.  e.  At  =  0),  In  order  to  perform  the  correlation,  the  analog  Doppler 
record  has  to  be  digitized  or  sampled  in  time  and  frequency.  For  this  purpose 
the.  time  and  the  frequency  axes  are  scaled  at  intervals  of  1/24  Hz  and  30  secs, 
respectively.  The  digitization  of  the  signal  is  accomplished  by  assigning  the 
Doppler  shifts  of  the  signal  either  J  or  0,  depending  on  whether  the  signal  is 
presort  or  not  in  the  given  cell  of  frequency-time  grid.  The  analog  and  digitized 
versions  of  the  Doppler  shift  records  at  the  3  stations  nave  been  shown  in  Figures 
2  and  7  respectively  (Sections  III  .  1  ind  IV),  The  course  digitization  might  have 
some  smearing  effect  on  the  corr  lation  functions  but  it  is  not  expected  to  cause 
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any  significant  difference  in  the  measurement  of  the  time  delays.  The  one-bit 
digitized  Doppler  values  are  denoted  as  Pj(t,  f),  I%(t,  f)  and  P^t,  f)  for  the 
3  stations  -  -  Ramey,  Arecibo  and  Mayaguez.  The  cross  correlation  function  for 
any  pair  of  signals  P^(t,  f)  and  Pg(t,  0  corresponding  to  a  time  lag  of  At  is 
expressed  as: 


P  A,  B<At> 


([PA<t,  f)  -  PA(t,  0]  [PB(t+At,  f)  -  PB(t,  f)]) 

2\>!  TJTr 


PA(t,f)  -  Pjt.f) 


1  1 


pA<'.  0  ■  V  <V*7> 


2\^l/2 


PB(t,f)  -  PR(t,  f) 


B' 


']>} 


I_f0  t=t2 

I  l 

f=  1  t=t. 


PA(t.  f) 


(-~‘o  *“*2+** 


PB(t’ f}=  f0  (t2-tj)  Z  Z  FB(t'f) 


f=  1  tst  |"r  At 


where  fg  is  the  total  number  of  frequency  scalings  and  (t^  -  tj)  is  the  length  of 
selected  window  for  correlation.  In  the  present  analysis  fg  =  24  and  (t^  -  tj)  = 

20  minutes.  The  c:oss  correlation  functions  p^At),  p^(At)  and  P23(At)  between 
the  pair*  of  stations  Ramay-Arecibo,  Ramey-Mayaguez  and  Arecibo-Mayaguez 
are  computed  for  At  varying  from  -10  to  +10  minutes  in  steps  of  one  minute.  In 
Figure  13  are  shown  the  plots  of  the  correlation  functions  for  six  different  cor¬ 
relation  windows  each  removed  from  the  other  by  5  minutes  over  the  total  test 
interval  14:08:30  -  14:53:30  U.  T.  The  plots  reveal  that  over  the  interval 
1408-1433  corresponding  to  the  first  two  correlation  windows  during  which  the 
major  complexity  of  the  Doppler  shift  occurred,  the  correlations  for  all  the 
3  pairs  of  stations  are  found  to  be  high.  During  the  later  periods,  however,  the 
correlation  is  found  to  remain  significant  only  between  Arecibo-Mayaguez  while 
for  the  other  two  pairs  the  correlation  has  degraded  rapidly  with  time.  Low 
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Figure  13.  Cross  Correlograms  for  Six  Successive  Positions  of  a 

20  Minute  Running  Window  in  the  Test  Interval  14:08:30  - 
14:53:30  U.  T.  on  June  12,  1967 


* 
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correlation  between  any  pair  of  the  stations  would  normally  mean  that  the  dis¬ 
turbance  front  is  small  in  dimension  compared  to  the  receiver  separation  or  that 
'he  shape  of  the  disturbance  is  changing  rapidly  compared  to  the  travel  time  of 
the  disturbance  from  one  receiver  to  the  other.  The  information  about  the  speed 
and  direction  of  the  disturbance  that  imposed  the  major  complex  signature  at 
1420  U.  T.  is  derived  f-om  the  first  correlation  window  which  furnishes  the 
correlation  coefficients  of  highest  significance.  The  velocity  is  calculated  from 
the  following  equations  using  the  station  geometry,  and  the  time  delays  derived 
from  the  correlation  analysis. 


-  cot  a 


■(£) 


cos  (a -9) 


whe  r  e : 

and  x  ^  are  the  distances  between  stations  1-2  and  1-3  respectively. 
T  and  are  corresponding  time  delays, 
a  is  the  angle  between  the  paths  joining  1-2  and  1-3. 

9  is  the  d’rection  of  propagation  of  the  disturbance  measured  counter¬ 
clockwise  from  the  line  joining  the  stations  1  and  3,  and 

S  i.  the  speed  of  propagation  of  the  disturbance.  The  velocity  and  the 
parameters  pertinent  to  the  calculation  are  listed  below. 


12 

1! 

km 

t  12  =  6.25  Min. 

S  =  112  m/s 

13 

=  32 

km 

Tj  ,  =  1.  50  Min. 

Bearing  =75°  west  of  north 
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B.  Ray  T racing 

1 .  Specification  of  the  Electron  Density  Distribution  in  the  Ionosphere 


The  observed  ionosphere  is  read  into  the  computer  as  a  two  dimensional 
(height-time)  array  of  plasma  frequencies.  The  interval  between  adjacent  points 
is  6  kilometers  in  height  and  33  seconds  in  time.  Once  in  the  computer,  the  time 
axis  of  the  array  is  converted  to  horizontal  distance  using  the  wave  speed  derived 
from  spaced  sounder  measurements,  the  plasma  frequencies  are  converted  to 
electron  densities  and  the  array  is  then  in  the  standard  form  for  ray  tracing. 

2.  Basic  Assumptions 

2.  1  Electron  Density  Between  Grid  Points  of  Input  Ionosphere 

The  electron  density  between  grid  points  is  found  by  linear  interpolation  in 
both  space  dimensions.  Figure  14  gives  a  physical  picture  of  the  interpolation  for 
a  typical  cell.  The  X-axis  represents  horizontal  distance,  the  Z-axis  represents 
altitude,  and  the  N-axis  represents  electron  density.  Linear  interpolation  amounts 
to  passing  a  plane  through  the  three  grid  points  in  the  X-Z-N  space. 

In  the  X-Z  plane  the  contours  of  constant  electron  density  are  straight  line® 
for  a  particular  cell.  Figure  15  shows  the  contours  within  two  adjacent  cells 
having  grid  point  densities  of  8,  9,  12  and  13  (arbitrary  units).  The  electron 
density  is  continuous  across  cell  boundaries.  The  magnitude  and  direction  of  the 
electron  density  grudien.  is  consta.it  within  each  cell  but  may  change  discontinu- 
ously  across  the  boundary. 

2.  2  Earth's  Curvature  and  Magnetic  Field 

The  earth's  curvature  ar.d  magnetic  field  are  ignored.  The  neglect  of 
curvature  is  of  no  practical  consequence  for  the  nearly  vertical  ray  paths  of 
interest.  The  neglect  of  the  magnetic  field  is  a  more  restrictive  assumption  but 
is  felt  to  be  justified  by  the  fact  the  time  displacement  between  disturbance 
signatures  observed  on  ordinary  and  extraordinary  rays  is  always  small  com¬ 
pared  to  the  wave  period  and  is  often  undetectable.  The  lateral  deviation  of  the 
rays  by  the  magnetic  field  is  therefore  of  secondary  importance.  Rays  traced 
using  the  true  operating  frequency  are  identified  as  ord’  ary  rays  for  the  purpose 
of  comparing  the  results  with  observation.  Extraordinary  rays  are  traced  using 
the  quasi-longitudinal  approximation  in  which  the  operating  frequency  is  reduces 
by  half  the  gyro-frequency. 
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Figure  15.  A  Representation  for  Typical  Electron  Density  Cells 


3.  Ray  Tracing  Method 

3.  1  Coo-dinate  Transformation 


Under  the  above  assumptions  rays  follow  parabolic  paths  within  each  cell.  ^  ^ 
The  parabola  differs  from  cell  to  cell  in  accordance  with  the  magnitude  and  direc¬ 
tion  of  the  electron  density  gradient,  the  operating  frequency,  and  the  initial 
orientation  of  the  ray.  The  equations  become  particularly  simple  when  the 
gradient  is  vertical,  the  zero  level  of  electron  density  is  at  the  ground,  and  the 
ray  is  launched  from  the  origin  of  the  coordinate  system.  To  take  advantage  of 
this,  the  ray  is  traced  thiough  each  cell  by  first  transforming  the  cell  from  its 
original  X-Z  coordinate  system  to  the  U-Y  .system  in  which  the  simple  equations 
are  valid.  The  transformation  consists  of  a  rotation  to  make  the  gradient  vertical, 
a  vertical  translation  to  make  the  zero  level  of  electron  density  lie  along  the 
horizontal  axis,  and  a  horizontal  translation  to  make  the  ray  pass  through  the 
origin.  Figure  16  illustrates  the  translation  for  a  typical  cell. 

3.  2  Ray  Tracing  Equations  in  the  Transformed  System 

Let  a  be  the  space  rate  of  change  of  electron  density,  normalized  to  the 
critical  electron  density  for  the  operating  frequency  used: 


In  MKS  units,  N  .  =  1.  24  x  10^  (f  ) ^  where  f  is  the  operating  frequency. 

crit  op  op  e.  / 

Similarly,  let  (3  be  the  time  rate  of  change  of  electron  density  normalized  to  N 

crit’ 


(3  -  (8N/8t)/Ncrit 


Then  the  ray  path  is  given  1  / 


sin 


2.  e. 


2  sin  6 


U  = 


0 


{cos^  -qV ^ 
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The  phase  path  length  from  the  origin  to  any  point  along  Le  ray  is 


P  =  [  ( 1+2  sin20Q)  cos  1+2  sin20Q  -  aV)  (cos -  a.^.}  2  ] 

The  instantaneous  Doppler  shift  imposed  along  the  ray  path  from  origin  to  ob¬ 
servation  point  is 

D  =  -(p/a>.)  [cos  ±  (cob^  -  aV)^^] 

where  X  is  the  operating  wavelength  in  free  space.  In  all  three  equations  the 
(+)  sign  is  taken  when  the  observation  point  lies  beyond  reflection  and  the  (-)  sign 
when  it  lies  before. 

3.  3  Procedure 

The  ray  is  traced  cell  by  cell.  The  first  step,  from  ground  to  ionosphere 
is  a  straight  line.  The  impact  point  on  the  bottom  of  the  ionosphere  determines 
the  first  cell  entered  by  the  ray.  The  cell  is  transformed  to  the  U-V  system  and 
the  ray  path  computed.  This  determines  the  point  at  which  the  ray  exits  from  the 
cell  and  its  orientation  at  that  point.  The  phase  path  length  and  the  Doppler  shift 
are  computed  from  the  origin  for  both  the  entry  and  exit  points  and  the  results 
subtracted  to  give  the  contribution  of  the  path  segment  within  the  ceil.  The  exit 
point  and  ray  orientation  as  the  ray  leaves  the  cell  constitute  the  entry  point  and 
ray  orientation  for  the  next  cell  and  the  process  is  repeated  until  the  ray  returns 
to  earth.  The  phase  path  length  and  Doppler  shift  contributions  from  each  cell 
are  added  as  the  ray  proceeds.  When  the  ray  impacts  at  the  ground  the  range 
from  the  transmitter  is  noted  and  the  take-off  angle  of  the  next  ray  modified 
accordingly.  After  several  trials  the  take-off  angie  required  for  the  ray  to  return 
to  the  transmitter  is  found.  The  ionosphere  is  then  moved  horizontally  with 
respect  to  the  transmitter  and  the  process  continued  until  a  Doppler  shift  versus 
time  recoid  has  been  synthesized  for  the  period  of  interest. 

4.  Cost 

The  calculations  involved  in  tracing  a  single  ray  require  about  one  second 
of  computer  time  (IBM  7044).  The  cost  per  ray  is  about  10  cer.fs. 
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C.  Wave  Interaction  With  Plasma 


An  atmospheric  wave  propagating  in  the  ionosphere  interacts  with  the 
ionospheric  plasma  and  causes  perturbations  in  the  electron  and  ion  densities. 

In  this  section  is  developed  a  simple  model  for  the  electron  density  perturbation 
in  response  to  an  acoustic  gravity  wave  propagating  in  a  non-isothermal  atmospher 
The  purpose  for  the  derivation  of  this  model  is  primarily  to  make  a  detailed  com¬ 
parison  of  the  theoretically  predicted  electron  density  fluctuations  with  those  ob¬ 
served  by  means  of  UHF  backscatter  measurements  described  in  an  earlier 
section.  The  interaction  between  the  neutral  gas  wave  motion  and  the  plasma 
is  governed  by  the  collisions  between  the  charged  and  the  neutral  particles  as 
well  as  the  geomagnetic  field  which  constrains  the  motion  of  the  charged  particles. 
In  fhe  region  above  about  150  km  both  the  ion-neutral  and  electron-neutral  col¬ 
lision  frequencies  are  much  less  than  the  respective  gyro-frequencies  of  the 
charged  particles.  Therefore,  for  motions  having  periods  in  excess  of  the  ion- 
neutral  collision  time,  the  only  effect  of  a  neutral  gas  motion  is  to  impart  to  the 
ionization  its  component  of  velocity  along  the  magnetic  field  line's.  Once  the 
electron  velocity  field  is  determined  from  the  neutral  gas  velocities  assuming  an 
appropriate  coupling  mechanism,  the  density  fluctuations  can  be  derived  by  solving 
the  electron  density  continuity  equation. 

1 .  Neutral  and  Electron  Gas  Velocities 

The  horizon^l  and  vertical  components  U  and  U  of  the  neutral  gas  per- 

X  z 

turbauv.  .«.y  are  expressed  in  complex  Fourier  form  from  wave  solutions 

of  the  linearized  equations  of  motion,  adiabatic  state  and  mass  conservation  Dy 
Hines^  ^  as: 

Ux  =  AX  exp  i(ut  -  Kxx  -K^z)  1(a) 

U  =  AZ  exp  i(wt  -K  x  -K  z)  1(b) 

Various  symbols  in  these  and  the  following  equations  in  this  section  carry 
the  same  meaning  as  given  by  Hines,  unless  otherwise  specified.  The  dispersion 
equation  relating  the  wave  numbers  and  the  frequencies  shows  that  when  Kx  is 
real,  say  K  =  k  ,  K  is  either  purely  imaginary  or  complex.  The  class  of  wavea 
which  are  called  internal  waves  which  exhibit  both  amplitude  and  phase  variations 
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in  vertical  direction  are  of  interest  in  studying  the  ionospheric  disturbances.  For 
the  case  where  the  viscous  damping  of  the  wave  amplitude  is  ignored,  the  dis¬ 
persion  relation  leads  to  K  =  k  +  and 

z  z  L  H 

z 

U  =  AX  exp  yvj-  exp  1  (ut-k  x  -  k  z)  2(a) 

X  ^ri  X  z 

=  AZ  exp  exp  i  (wt-k^x  -  k^z)  2(b) 

where  A  is  a  constant  and  X  and  Z  are  the  polarization  factors  given  as 


X  =  oo  k  C2  [k  -  i  (1  -  v/2)  g/C2]  3(a) 

X  z 

Z  =  W3  -  w  k  2  c2  3(b) 

x 

The  wave  number  k  , which  is  real  in  equations  (2)  in  the  absence  of  viscous 

z  i 

dampingjWOuld  become  complex  (k?  =  k*)  when  the  viscous  effects  are  taken  into 

consideration.  Pitteway  and  Hines^^  obtained  a  dispersion  relation  considering 

viscous  damping  which  may  be  used  to  determine  the  complex  k*.  In  the  low  fre- 

quency  case  such  as  of  interest  in  the  present  analysis,  some  valid  approximations 

in  the  dispersion  relation  yield  the  folic  wing  relation  for  k*: 

z 


k 


z 


+  k 


2  2  2 
2  k  «//</ 
*  g 


1  +  Vl-iD 


(4) 


where : 

D  =  4t|  k  2  U  2/w3  (5) 

x  g 

D  would  have  to  be  taken  as  2.  26  times  that  given  in  equation  (5)  if,  as  necessary, 
the  effect  of  thermal  conduction  is  taken  into  account  (Hines,  Personal 
Communication). 
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On  separating  the  real  and  imaginary  parts  in  Equation  (4)  we  have 


r  7 ,  2  2 

2  k  w 

...  x  g 

1  I - 2  \'/2  1 

VI  +  D  -  1  2  1 

2 

L  w  D 

i  2  /  *  «2J 

Denoting  the  real  and  imaginary  parts  in  the  above  equation  as  'a'  and  'b' 
respectively  we  have  for 


k 

z 


i. 


± 


|^aZ+b2+a 


+  i 


l/2_ 

2  j 

(6) 


The  choice  of  sign  for  kz  is  governed  by  the  requirement  that  the  amplitude  of 
the  wave  is  attenuated  progressively  in  the  positive  z  direction  due  to  viscous 
damping.  The  complex  k*  which  may  be  expressed  as  'k  -  iy)  where  k  and  x 
are  real  modifies  the  perturbation  velocity  components  U  and  U  as: 


U 

x 


AX  exp  z  (-j^-  -  X)  exp  i  (wt  -  kxx  -  kzz) 


7(a) 


U 

z 


-  AZ  exp  z  ( 


1 

2H 


-  x)  exp  i  (wt  -  k  x  -  k  z) 


7(b) 


The  polarization  factor  X  is  now  modified  in  that  the  kz  in  equation  3(a)  is 
replaced  by  kz>  Z,  however,  remains  the  same  as  in  3(b).  The  component  of 
the  neutral  gas  velocity  along  the  geomagnetic  field  which  is  effective  in  driving 
the  plasma  can  be  expressed  for  the  velocity  and  magnetic  field  configuration 
shown  in  Figure  17  as: 
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Figure  17.  The  Neutral  Gas  Velocity  and  the  Geomagnetic  Field 
Configuration 

V-.  =  (U  cos  0  cos  I  -  TJ  sin  I) 
n>  x  z 

where  9  is  the  angle  of  direction  of  U  measured  west  of  north  and  I  is  the 
magnetic  dip  angle. 

The  horizontal  and  vertical  components  V  and  V  of  the  electron  gas  in 

x  z 

the  directions  of  U  and  U  are  respectively: 

x  z 

V  =  V„  cos  I  cos  9 
x  B 

Vz  =  -VB  sin  I 

Substituting  for  Vg  and  then  expanding  for  and  we  get: 
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2  2 

X  cos  9  cos u  I  -  Z  sin  I  cos  I  cos  9 


A.  exp  z  (  2jT  -  X  ) 


exp  i  (ut  -  k^x  -  k^z) 


8(a) 


Z  sin^  I  -  X  sin  I  cos  I  cos  9 


a  /  1  \ 

A  exP  1  ( liT  "  x  ) 


exp  i  (wt  -  k  x  -  k  z) 


8(b) 


on  separating  the  rial  and  imaginary  parts  in  the  amplitude  factors  the  equations 
8(a)  and  8(b)  may  be  written  in  the  following  form: 


Vx  =  (P+iQ)  A  exp  z (~-  -  x)  exp  i  (ut  -  kxx  -  kzz) 


9(a) 


Vz  =  (R+iS)  A  exp  z(-=p^-  -  x)  exp  i  (wt-k^x  -  kzz) 


9/1  ' 


where 

P  =  wk  k  G2  cos2  6  cos2  I  -  (u^  -  wk^  C2)  sin  1  cos  I  cos  9 

X  Z  X 

Q  =  -uk^  t  [X  C2  +  (1  -y/2)  g]  cos^  9  cos^  I 
R  =  (w^  -  w  k  C2)  sin  I  -  wk  k  C2  sin  I  cos  I  cos  9 

X  X  z 

and 


S  =  wk  .  tx  C2  +  (1  -  y/2)  g J  sin  I  cos  I  cos  9 


The  velocities  V  and  V  of  the  electron  is  can  now  be  used  to  calculate  the 
x  z 

electron  density  perturbations  from  the  continuity  equation  for  the  electron  density. 
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2 .  Continuity  Equation 


The  continuity  equation  for  both  electrons  and  ions  in  the  ionosphere  where 
their  densities  arc  assumed  to  be  equal  may  be  expressed  m  the  form: 


$  +  V  •  [N  (Vd  +  Ve)J  =  Q  -  L  (10) 

where  Q  and  L,  are  the  production  and  loss  rates 

N  =  electron  (or  ion)  density 

=  diffusion  velocity  of  the  ionization 
Vg  =  drift  velocity  due  to  electric  fields 

L  has  the  form  of  f3N  in  the  F  region  and  ciN^  in  the  E  region. 

In  the  present  analysis  L  is  represented  by  (3N  by  adopting  an  eqi  ivalent 
for  P  for  the  heights  below  about  180  km  where  the  loss  is  no*  of  purely  attachment- 
like. 

The  electron  density  14  is  perturbed  under  the  i  ..  aence  of  an  acoustic- 
gravity  wave  so  that  the  perturbed  density  N  may  be  written  as: 

.  i(ut-k  x-k  z) 

N  (z)  =  N(z)  r  n(z)  =  N(z)  +  n^r.)  e  x  z  (11) 

The  continuity  equation  for  the  perturbed  ionosphere  may  be  written  as: 


9N 


9t 


+  V 


[NJ(V  +  V,  +  V  )]  =  Q  -  pN1 


(12) 


whv're  V  is  the  wave  induced  velocity  of  the  ionization.  Subtracting  equation 
{ lu)  from  (12)  we  have  the  following  equation  fo*-  density  perturbation  n. 


+  V  •  t  N1  V  +  niV  ,  +  V  )  j  =  -Pn 
at  a  e 


(13) 
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Assuming  that  the  density  perturbation  n  is  much  smaller  than  the  ambient 
electron  density  N  and  also  the  contribution  due  to  diffusion  and  the  electric  fields 
is  unimportant  for  the  wave  periods  of  interest/^’  ^  the  above  equation  is 
simplified  to: 


+  pn  =  -V-  (NV)  (14) 

Neglecting  the  horizontal  gradients  in  the  ambient  electron  density  and  expanding 
various  terms  of  equation  (14)  we  have. 


n(p+i  u:) 

r  av 

M  X  1 

=  - N  Lax  4 

av 

av 

X 

ilc  V  1  Z 

dx 

xx’  3z 

"  -  v  4^ 

z  dz 


-  ik  J  V 
z  /  z 


n  [p+iUJ  -  Vx[lNkJ+Vi  [n,'v  --ji-tlk,)  -!£-] 


=  A  exp  2  (- 


X)  exp  i  (wt  -  kxx  -  kzz) 


where 


j(P+  iQ)  iNkx+  (R+iS)  N(X  -  -gjj-  +  ikB)  -  j 

A  exp  z  (-5—  -  x)  exP  i(«Jt  -  k  x  -  k,  z)  [f  +  i*] 

fail  X  Zt 


1  -  N  [Rx  -  ifr  - Q  kx  -  Skz 


]  -«(-&) 


(  =  N  Pk  +  Rk  +  Sx  -  TET 

lx  z  ZH 


-  s  (If) 
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A  exp  z  ( 


I 

2H 


-  x)  e*p  i(wt-k  x-k^z) 


[  (fp  +  geo)  +  i(  ip  -fu>)] 


(15) 


The  above  equation  can  be  used  to  compute  the  electron  density  perturbation 
n  as  function  of  height  and  time  at  any  given  location  knowing  the  following 
information: 

1.  The  wave  parameters  ~  and  (k^  may  be  computed) 

2.  A  model  lor  the  ambient  electron  density  N(z,  t) 

3  An  atmospheric  model  to  calculate  H,  x  and  various  other  parameters 
needed  for  f  and  g. 

4.  A  model  for  the  loss  coefficient  (3(z),  and 

5.  The  constant  A. 

A  computer  program  to  calculate  the  electron  density  perturbation  n(z,  t) 
is  being  de  eloped  with  the  observed  HF  and  UHF  data  providing  the  wave 
parameters  and  the  ambient  electron  density  model.  An  atmospheric  model 
appropriate  t0  the  period  of  observation  is  selected  from  CIRA  (1965).  For  the 
attachment  coefficient  p(z),  a  model  originally  suggested  by  Ratcliffe  et  al  ’ 
and  used  extensively  in  later  studies  will  be  adopted  here.  The  constant  A  will 
be  determined  by  stating  the  initial  value  for  the  perturbation  which  would  be 
provided  from  the  observations.  It  is  now  required  that  the  density  perturbation 
be  extracted  from  the  observations  for  a  comparison  with  the  theory.  The  in¬ 
coherent  backscatter  measurements  furnish  the  wave  perturbed  electron  density 
profiles  n\z,  t).  The  wave  perturbation  in  the  observations  is  averaged  out  by 
performing  the  running  averages  of  N*(z,  t)  with  a  Lime  scale  equal  to  the  pre¬ 
dominant  wave  period  which  leads  to  the  ambient  electron  density  profiles  N(z,  t). 
The  difference  of  N  (z,  t)  and  N(z,  t)  gives  the  perturbation  profiles  n(z,  t). 
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